Massless Leggett Mode in Three-band Superconductors with Time-Reversal-Symmetry Breaking 
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The Leggett mode associated with out-of-phase oscillations of superconducting phase in multi-band super- 
conductors usually is heavy due to the interband coupling, which makes its excitation and detection difficult. We 
report on the existence of a massless Leggett mode in three-band superconductors with time-reversal-symmetry- 
breaking (TRSB). The mass of the Leggett mode is small close to the TRSB transition and vanishes at the 
transition point, and thus locates within the smallest superconducting energy gap, which makes it stable and 
detectable. The mass of the Leggett mode can be measured by Raman spectroscopy. The thermodynamic 
consequences of this massless mode and possible realization in iron-based superconductors are also discussed. 
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Introduction - Spontaneous breaking of a continuous symme- 
try and the associated low-energy collective excitation gov- 
ern the physical properties in many systems ranging from 
condensed matter physics to particle physics. Supercon- 
ductivity emerging as the spontaneous breaking of the t/(l) 
gauge symmetry supports a massless excitation known as 
Bogoliubov-Anderson-Goldstone (BAG) bosonfl, 2]. Cou- 
pled with electromagnetic field, the BAG boson becomes the 
massive plasma mode due to the Anderson-Higgs mechanism. 

Because of the discoveries of MgB2l3)] and iron 
pnictides|4], it is now accepted that multi-component super- 
conductors are ubiquitous. Multi-band superconductors are 
not straightforward extensions of the single-band counterpart, 
novel features arise insteadpM.KJO. A famous example is the 
Leggett mode (LM) in two-band superconductors associated 
with the collective oscillation of superconducting condensates 
between different bands, as schematically depicted in Fig. [T] 
with the mass proportional to the interband coupling. [51] In 
2007, Blumberg et. at. reported the observation of the LM in 
MgB2 with the Raman spectroscopy [9]. The mass of the mode 
lies between the two superconducting energy gaps, consistent 
with the theoretical calculations! 1 II- The LM in MgB2 there- 
fore decays into quasiparticle continuum associated with the 
band of smaller energy gap. The heavy LM in MgB2 has also 
been observed in point-contact transport measurements |1 2]. 

For iron-based superconductors, many studies have re- 
vealed the sign-reversal pairing symmetry between different 
bands il3l4l6ll . The system of more than three bands is some- 
how frustrated, and under appropriate conditions there may 
exist time-reversal-symmetry-breaking (TRSB) states even 
with conventional s-wave pairing symmetry, which involve 
non-trivial phase differences (i.e. 6(p ^^ or n) among su- 
perconducting gaps|6, 17-19]. With the new TRSB transi- 
tion below Tc, the transition temperature of superconductiv- 
ity, the spectrum of collective excitations and thus low-energy 
physical properties of the superconductors should be modified 
significantly. It was reported that the LM may exist below 
the two-particle continuum in iron-based superconductors un- 
der appropriate conditions, 12011 and that the mass of the LM 
may be reduced in some dynamical classes of multiple inter- 
band Josephson coupling in three-band superconductors. 112 1 1 



FIG. 1 . Frustrated interband scatterings force Cooper pairs in differ- 
ent bands carry different phases, which results in interband Joseph- 
son currents. Two dynamical modes associated with superconduct- 
ing phases in three-band superconductors: the LM, where one of the 
three phases stays unchanged while the other two oscillate out-of- 
phase, becomes massless at the TRSB transition (left), and the BAG 
mode, where all the three phases rotate in the same direction during 
the propagation of plasma wave in space (right). 



We note that in another TRSB superconducting systems with 
mixed-symmetry order parameters with nodes such as li -i- is, 
a massive LM in the TRSB state was found in Ref. [22[]. With 
all these recent progresses in mind, we ask a question which 
is of fundamental interest: what are the effects of frustration 
on the LM and is it possible to have a massless LM mode? 

In the present work, we demonstrate that the mass of the 
LM can be reduced significantly and even vanishes at TRSB 
transition upon sweeping interband coupling or density of 
states in multi-band superconductors. It is shown that the LM 
of vanishing mass can be detected by Raman scattering, which 
also serves as a smoking gun evidence for the TRSB transi- 
tion. The appearance of massless excitation modifies super- 
conductivity properties qualitatively, such as the power-law 
dependence of the specific heat (SH) on temperature instead 
of the conventionally exponential one for full-gapped systems. 
It is found that several recent experiments on iron-based su- 



perconductors can be explained by the existence of massless 
LM. 

Leggett mode - The Hamiltonian for three separate pieces of 
the isotropic Fermi surface can be written as 



jJ 



(1) 



where tfrj^ (if'io-) is the electron creation (annihilation) opera- 
tor in the /-th band with the dispersion e/(k) and the chemical 
potential fi and spin index cr. Vji is the intraband for I = j and 
interband for / ^ j scattering respectively, which can be either 
repulsive or attractive depending, for instance, on the strength 
of the Coulomb and electron-phonon interaction. The inter- 
band repulsion may cause frustration of the superconductivity 
in different bands and results in TRSBdallTI]. Introducing the 
Nambu spinor operator ^'j - (tffji, lA i )^ ^"d the energy gap 
Aj through the Hubbard-Stratonovich transform, we arrive at 
the following action in the imaginary time representation after 
integrating out the fermionic fields [i23ll 



(2) 



JdTd'rf.^jgjiA'i-Yj'^ringj' 



jJ 



with f = y ' and the Gor'kov green function 



^J' = - 



dr + (ej-fi) -Aj 



The superconducting energy gaps at T = are given by 

3 



Y^A,gij = Nj{0)Ajsmh- 



1=1 



hcUcj 



(4) 



with A^/(0) the density of states (DOS) at the Fermi surface in 
normal state. Here a>cj is a cutoff frequency and depends on 
the pairing mechanism. For electron-phonon coupling, cJcj is 
the Debye frequency. 

For demonstration of our basic idea, we take a set of sim- 
plified interband couplings ll24ll 
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and assume that the DOS A^ and a>c are identical for three 
bands. [17] The massless LM, however, is not restricted to the 
specific choice of g as discussed later. Here gij > corre- 
sponds to a repulsive interaction. We take Ai as positive real, 
and A2 = Ae'^, A3 = Ae^"^ because they are symmetric under 
the condition of Eq. (|5]l. Hereafter we take ftoJc as the unit for 
A,. 

For a small 77, the interband repulsion ^12 and go dominates 
and the system takes if - n. For a large 77, a state with finite 
phase diff'erence between A2 and A3 appears, corresponding to 
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FIG. 2. (color online). Amplitudes and phases of order parameters at 
TRSB phase transition, in (a) and (b) as a function of 77, and in (c, d) 
as a function of DOS N\ V of the first component. Ai is taken as real 
and positive. In (a) and (b), an identical DOS NY = 0.5 is taken for 
the three bands and a = 2 in Eq. (5). In (c) and (d), N2V = 0.5 and 
N3V = 0.4, tt = 2 and 77 = I (see also Eq. (17)). In the TRSB regime, 
there are two degenerate ground states (Ai,A2, A3) (solid lines) and 

refer to the same state without TRSB. 



A*, A*) (dashed lines). The two solid lines for NiV > 0.64 in (d) 



astateof TRSB where (Ai,A2, A3) ^ (Ai,A2, A3)*, even apart 
the common phase factor. In the TRSB state the energy gaps 
are given by 



'''"^(^;^)-'"^* = "-''(¥^) <« 



and cos ip - -Ai/drjA). The system undergoes a second- 
order TRSB transition at 77^. given by 77c = Ai(77c)/[2A(?7c)], as 
shown in Fig. |2](a) and (b). 

We proceed to investigate phase fluctuations at the TRSB 
point, where the amplitudes of the superconducting gap can 
be considered as rigid. For this purpose, we perform the fol- 
lowing gauge transformation which separates the phase and 
amphtude of gap illll25r 



Aj -^ lAjle'"! and *P/t, r) 



'0 ,-«,/2 l^/T,'-). (7) 



and derive the action for the phase fluctuation 
5 = r dTd'rY, |A/lgo|A,k'(«'-»'> - 2 Ti- [in {&]' - ^j)] 



(8) 



where S,- = -^i{V^ej + N0jV)cro + [i'-^ + £-(V0,) ]cr, 
with (Tj being the Pauli matrices, ctq the unit matrix and nij 



the electron mass II26L 12711 . From this action, one can obtain 



the time-dependent nonlinear Schrodinger Lagrangian for the 
phase fluctuations i28i 12911 . Considering small phase fluctua- 
tions around the saddle point (pj = 6j - ipj and expanding S up 



to the second order in ^j, we have||24 1 

S^ h] -\Yj jd'q^i-^i, -qVM4>(a,, q) (9) 



with ^(Q;, q) = [0i(Q;, q), (piiQ'i, q), hi^h q)Y and 
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with Di - 8AiAcos^/y and D2 - 877A cos(2^)/y with (p = 
ip2-f\ - (fii -tp3- ^1 - llnksT and the excitations are bosons. 
In the hydrodynamic Hmit at T = 0, the dissipation is absent 
and Pj = INi-Qr' + l/3v^.q^) after the analytical continuation 
/Q/ <— fi + iO^. From DetM - 0, we obtain the dispersion 
relations 
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The first mode is the massless BAG mode corresponding 
to the uniform rotation of phases. The second and third are 
the LM Ql- and Ql+ in the present three-band system. Espe- 
cially, the mode Ql corresponds to the dynamics of the rel- 
ative phase ip23 between the gaps of A2 and A3, and becomes 
massless at the TRSB transition depicted in Fig. [3] One may 
regard if 23 as the order parameter for the TRSB transition. As 
it increases continuously from at the transition, the associ- 
ated fluctuations become massless at the TRSB transition. 

The magnetic field can be introduced into S^ through the 
standard replacement V0/ -^ V(f)i - 2;rA/<l>o with Oq the flux 
quantum and A the vector potential. In this case, it is more 
convenient to describe the phase fluctuations in terms of (pi, 
012 = 01-02 and 013 =01-03. 01 describes the BAG 
mode, and 0i2 and 0i3 correspond to the LMs. The gauge 
field couples with 0i in the form (V0i - 27rA/<l'o)- One may 
integrate out 0i, resulting in the massive plasma mode due 
to the Anderson-Higgs mechanism. In contrast to the BAG 
mode, the LMs remain massless at the TRSB transition since 
012 and 013 are decoupled from the gauge field A. 

In stark contrast to conventional symmetry-broken systems, 
there exist the stable LMs both before and after TRSB transi- 
tion, because the relative phase between different condensates 
is fixed in both the states with and without TRSB. 

Raman scattering- Interband scatterings do not involve the 
gauge field, thus the LMs do not respond to a magnetic field. 
However, the LMs are coupled indirectly with the electric 
field through the charge density, which renders it detectable 
by the Raman spectroscopy through the inelastic scattering of 
photon with the charge density ||30. - .33.1 . The interaction be- 
tween the incident photon and the charge can be modeled as 

piT, ?) = Z Z JjWilf]^ (r, k+i) ^ja (r, k - f ), where yj{k) 

j=\ k.cr 




FIG. 3. (color online). Dependence of the masses of LMs on the 
interband coupling rj. Here NV = 0.5 and a = 2, and the masses are 
in units of w^. 



is the scattering coefficient determined by the polarization 
of the incident and scattered photon. In the following, we 
derive the experimentally measurable Raman response func- 
tion ;^fpp(T - r', q) - - {TrPir, q)p{_T' , -q)) with T^ being time- 
ordering operator. We introduce a source term coupled with 
p, Hj{t) = - 2 p{t, q)J{T, -q) because xpp can be computed 

by the linear response theory with respect to J. The effective 
action in the presence of incident photon reads ll24|] 



S^ ^dTd\Yj^,g,jA)-Y,'^v\n{Q-j}+g]') (14) 

with ^J J = -yi{k)J{T,-q)o-3. For a weak incident wave, we 
may neglect the fluctuations of the amplitude of the order pa- 
rameters, and the fluctuations for the superconducting phase 
acquires a form S = Sa + Sj, with S^ defined in Eq. ^ and 



Sl = ^nJ(q)Zj(q)<p]i-q) 
+J(-q)Zj(-q)<Pj(q) + J(q)J(-q)n'j^^], 



(15) 



where Zj(q) - Aj[-sm(pjYV.j^{q) - cos ^^H'^ 32(^)1 '^^'^ 
^A4) = A^[-sin^jn^j3(g') - CQ^ifjIlT^^iq)]. The 
polarization functions are defined as 11'^'^, ,11^ , = 

1 /(L'/3) Zn j d'kYj,,,,, [yj{k + l)yj(k - f ), yj{k + f )] . 

Integrating out the fluctuations 0y, we then obtain the cor- 
relation function 



,p(/a q-0)-Yu K33 - Z;[M-'],vZj} . (16) 



Xpp 



The first term gives the resonant scattering at Q = 2Aj and 
the second term accounts for the resonance with the LMs, 
as depicted in Fig. |4] When the energy shift of the photon 
matches the energy of the LMs, M"' becomes singular and 
gives 6 peaks in the spectroscopy. In reality, the delta-function 
peaks are rounded by both damping effect and interactions be- 
tween Leggett bosons when the oscillations of the LM become 




strong frustration interactions among different bands. Here we 
wish to observe that the TRSB transition can be induced not 
only by interband coupUng but also by DOS Nj{Q). In order 
to demonstrate this we derive the TRSB solution to Eq. (HI 
under a general coupling matrix g. Complex gap functions as 
solution to Eq. (HJi appear when there is only one independent 
vector in the matrix g - g' , with g'.. - Nj{0) smh~^(hcori/\Ai\) 



Raman shift 



and otherwise. From this constraint we obtain ll24ll 



FIG. 4. (color online). Schematic view of the Raman response 
in three-band superconductors with TRSB. The finite line- width of 
peaks is due to damping and interaction between Leggett bosons. 
The background at energy larger than 2Ai is due to the quasiparticle 
excitations. 



strong, which are absent in Eq. ( fTSI l. Although the response of 
a genuinely massless LM is hidden into the elastic scatterings, 
it can be traced out clearly if one changes ;; systematically and 
generates LM of small mass, which can be achieved by car- 
rier doping because the interband scattering is renormaUzed 
by the DOS as in Eq. dUi. 

Discussions- At T > 0, the Landau damping by quasipar- 
ticles sets in and the lifetime of the LM decreases. A local 
time-dependent equation for the phase fluctuations does not 
exist due to the singularity of the DOS in the superconduct- 
ing state[34]. In the vicinity of T^, the dynamics of super- 
conductivity can be described by the standard time-dependent 
Ginzburg-Landau equation ll34ll . In this region, the lifetime 
of the LM is much smaller than the inverse of its energy 
due to the severe damping by quasiparticles, therefore there 
is no well-defined Leggett excitations. Nevertheless in the 
static case, the massless feature manifests as the divergence 
of the characteristic length for the relative phase variation in 
the vicinity of the TRSB 1 19]. 

Let us discuss the applicability of our results to the iron- 
based superconductors. In order to demonstrate the mass re- 
duction of the LM by TRSB, we adopt a simple and general 
BCS-like Hamiltonian in Eq. ([TJ. An implicit expectation be- 
hind this treatment is that more realistic models would merely 
lead to quantitative corrections. It seems that this simplifica- 
tion is not far from the situation in some iron-based super- 
conductors, since the .s-wave with sign-reversal (s+) pairing 
symmetry is favored by many experiments 11 1 5[ 1 1 611 . It also be- 
came clear recently that the interband hopping in single parti- 
cle channels using a more realistic tight-binding model gives 
additional contribution to the interband Josephson coupling, 
and that the Hunds interaction only gives higher order correc- 



hojcj sinh[(g jjgki - 8jk8ji)/NjiO)gki] ' 



(17) 



tion to the LM 1,20.1 . It was also shownl,35,l that the sipairing 
can result from the moderate electronic correlations l36l - l40ll in 



iron-based superconductors, thus electronic correlations prob- 
ably do not hamper much the massless Leggett mode. 

The reason that no direct experimental observation on the 
TRSB state in iron-pnictide superconductors has been re- 
ported to date may come from its requirement on sufficiently 



with j + k i^ l.li is easy to see that to find further the phases of 
the gap functions is equivalent to forming a triangle with the 
three segments \Aj\lgku which is possible when and only when 
\Aj\lgki + \Ak\lgji > |A;|/^j<: for all the three combinations. The 
phase transition from a TRSB state to a state without TRSB 
takes place when one of the above inequalities is broken, for 
example |Ai|/g23 = lAaJ/gn + lA^l/gn- The results for DOS- 
driven TRSB transition are displayed in Fig. |2](c) and (d). 
There are two TRSB transitions and the TRSB state is realized 
in a finite region of DOS. Therefore, experimentally one can 
tune Nj(Q) by careful chemical doping, which hopefully will 
drive the system to the TRSB transition. 

Although the massless LM does not change magnetic prop- 
erties of the system, it results in qualitatively different thermo- 
dynamic behaviors of .s-wave superconductivity. For the SH, 
the contribution due to quasiparticles atT -^ Tc depends ex- 
ponentially on temperature (A/fe^T^)^^^ sfi-pi-A/ksT) for fully 
gapped superconductors. The contribution of the massless 
Leggett excitations can be obtained analytically by treating 
the Leggett bosons as free quantum gas. The contribution is 
of power-law temperature dependence T^, which can be de- 
tected experimentally. 

It is worth noting that a T^ dependence of the SH in iron- 
base superconductors after subtracting the residue electronic 
contribution (linear in T) and phonon contribution (also T^ 
dependence) has been reported in several experiments, |4lj- 
43]; fully gapped order parameters are inferred from measure- 
ments for the dependence of electronic SH on magnetic field, 
which excluded the possibility of gap function of line node. 
Actually, in Ref. 1.4Z1 . the authors suggested that the addi- 
tional T^ contribution might be due to some bosonic modes. 
These experimental observations can be naturally explained 
by the existence massless LM. Additional measurements such 
as the Raman spectroscopy on similar samples ll4ll - l43ll are 
much anticipated which may well be in the vicinity of the 
TRSB transition. 
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